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Abstract 

Oh ' We consider the inviscid limit for the two-dimensional incompressible Navier-Stokes equa- 

i tion in the particular case where the initial flow is a finite collection of point vortices. We 

suppose that the initial positions and the circulations of the vortices do not depend on the 
viscosity parameter i/, and we choose a time T > such that the Helmholtz-Kirchhoff point 
vortex system is well-posed on the interval [0,T]. Under these assumptions, we prove that 
the solution of the Navier-Stokes equation converges, as — > 0, to a superposition of Lamb- 
Oseen vortices whose centers evolve according to a viscous regularization of the point vortex 
system. Convergence holds uniformly in time, in a strong topology which allows to give an 
^ ' accurate description of the asymptotic profile of each individual vortex. In particular, we 

, compute to leading order the deformations of the vortices due to mutual interactions. This 

' allows to estimate the self-interactions, which play an important role in the convergence 

^ . proof. 

od 

^ ■ 1 Introduction 

o ■ 



It is a well established fact that coherent structures play a crucial role in the dynamics of 
two-dimensional turbulent flows. Experimental observations [lOj and numerical simulations of 
decaying turbulence |301 l31j reveal that, in a two-dimensional flow with sufficiently high Reynolds 
number, isolated regions of concentrated vorticity appear after a short transient period, and 
persist over a very long time scale. These structures are nearly axisymmetric and behave like 
point vortices as long as they remain widely separated, but when two of them come sufficiently 
close to each other they get significantly deformed under the the strain of the velocity field, and 
the interaction may even cause both vortices to merge into a single, larger structure |29i I42j . It 
thus appears that the long-time behavior of two-dimensional decaying turbulence is essentially 
governed by a few basic mechanisms, such as vortex interaction and, especially, vortex merging. 

Although these phenomena are relatively well understood from a qualitative point of view, 
they remain largely beyond the scope of rigorous analysis. Vortex merging, in particular, is 
a genuinely nonperturbative process which seems extremely hard to describe mathematically, 
although it is certainly the key mechanism which explains the coarsening of vorticity structures 
in two-dimensional flows, in agreement with the inverse energy cascade. The situation is better 
for vortex interactions, which may be rigorously studied in the asymptotic regime where the 
distance between vortices is much larger than the typical core size, but complex phenomena 
can occur even in that case. Indeed, numerical calculations ^29j and nonrigorous asymptotic 
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expansions [50^ I49j indicate that vortex interaction begins with a fast relaxation process, during 
which each vortex adapts its shape to the velocity field generated by the other vortices. This first 
step depends on the details of the initial data, and is characterized by temporal oscillations of 
the vortex cores which disappear on a non- viscous time scale. In a second step, the vortices relax 
to a Gaussian-like profile at a diffusive rate, and the system reaches a "metastable state" which 
is independent of the initial data, and will persist until two vortices get sufficiently close to start 
a merging process. In this metastable regime, the vortex centers move in the plane according 
to the Helmholtz-Kirchhoff dynamics, and the vortex profiles are uniquely determined, up to a 
scaling factor, by the relative positions of the centers. 

Prom a mathematical point of view, a natural approach to study vortex interactions is to 
start with point vortices as initial data. After solving the Navier-Stokes equations, we obtain in 
this way a family of interacting vortices which, by construction, is directly in the metastable state 
that we have just described. In particular, as it will be proved below, we do not observe here 
the oscillatory and diffusive transient steps which take place in the general case. Point vortices 
can therefore be considered as well prepared initial data for the vortex interaction problem. 

With this motivation in mind, we study in the present paper what we call the viscous 
N -vortex solution^ namely the solution of the two-dimensional Navier-Stokes equations in the 
particular case where the initial vorticity is a superposition of point vortices. Por a given 
value of the viscosity parameter this solution is entirely determined by the initial positions 
xi, . . . , xat and the circulations ai, . . . , oat of the vortices. It describes a family of interacting 
vortices of size which are therefore widely separated if v is sufficiently small. Our 

main goal is to obtain a rigorous asymptotic expansion of the A-vortex solution in the vanishing 
viscosity limit, assuming that vortex collisions do not occur. As was already explained, this 
problem is physically relevant, but it also has its own mathematical interest. Indeed, it is known 
that the two-dimensional Navier-Stokes equations have a unique solution, for any value of 
when the initial vorticity is a finite measure [TSJ, but computing the inviscid limit of rough 
solutions is a very difficult task in general, due to the underlying instabilites of the Euler flow. 
Surprisingly enough, although point vortices are perhaps the most singular initial data that can 
be allowed for the Navier-Stokes equations, the inviscid limit appears to be tractable for the 
A- vortex solution, and provides a new rigorous derivation of the Helmholtz-Kirchhoff dynamics 
as well as a mathematical description of the metastable regime for interacting vortices. 

In the rest of this introductory section, we recall a global well-posedness result for the 
two-dimensional Navier-Stokes equations which is adapted to our purposes, we introduce the 
Lamb-Oseen vortices which will play a crucial role in our analysis, and we briefly mention the 
difficulties related to the inviscid limit of rough solutions. Our main results concerning the 
A^- vortex solution will be stated in Section [21 and proved in the subsequent sections. 

The incompressible Navier-Stokes equations in the plane have the following form: 

— + {u-\/)u = i^Au-Vp , divn = 0, (1.1) 

where m(x, t) G denotes the velocity of the fluid at point x S and time t > 0, and p{x, t) S M 
is the pressure inside the fluid. The only physical parameter in (jl.ip is the kinematic viscosity 
> 0, which will play an important role in this work. Por our purposes it will be convenient 
to consider the vorticity field u){x,t) = diU2{x,t) — d2Ui{x,t), which evolves according to the 
remarkably simple equation 

^ + (li • V)lo = i^Auj . (1.2) 
Under mild assumptions, which will always be satisfied below, the velocity field u{x, t) can be 
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reconstructed from the vorticity t) via the two-dimensional Biot-Savart law: 



u{x,t) = — ^-c{y,t)dy, xGR2, (1.3) 



2vr7R2 \x - 

where, for any x = {xi,X2) € M^, we denote X"*- = {—X2,xi) and |xp = xf + X2- 

Let A^(R2) be the space of all real- valued finite measures on R , equipped with the total 
variation norm 

||/i||tv = sup{(/u„, Co(R2) , UWl^ < 1} . 

Here (/>) = /]g2</'d^, and Co(R^) denotes the space of all continuous functions (/> : R^ — > R 
which vanish at infinity. We say that a sequence {/in} in AI(R^) converges weakly to fi £ Ai{M?) 
if {nn, (t>) — *■ (a*) as n ^ oo for all (f) G Co(R^). Weak convergence is denoted by /i„ ^ /i. 

Our starting point is the following result, which shows that the initial value problem for 
Eq. (jl.2p is globally well-posed in the space M. ( 



Theorem 1.1 ^15j Fix v > Q. For any initial measure fi G A4(R^), Eq. (11. 2p has a unique 
global solution 

a; G C°((0,oo),L^(R2)nL°°(R2)) (1.4) 
such that ||lij(-, t)||^i < ||/i||tv for all t > 0, and uj{-,t) ^ /u as t ^ 0+. 

Here and in what follows, it is understood that u is a mild solution of (jl.2p . i.e. a solution 
of the associated integral equation 

uj{t) = e^'^fi - div(e^(*-^)^u(s)u;(s)) ds , t > , (1.5) 

where e*^ denotes the heat semigroup. The existence of a global solution to (jl.2p for all initial 
data in 7W(R-^) has been established more than 20 years ago by G.-H. Cottet [9j, and inde- 
pendently by Y. Giga, T. Miyakawa and H. Osada [21]. In the same spirit, the later work by 
T. Kato [27] should also be mentioned. In addition to existence, it was shown in [211 [27] that 
the solution of (II. 2p is unique if the atomic part of the initial measure is small compared to the 
viscosity. This smallness condition turns out to be necessary if one wants to obtain uniqueness 
by a standard application of Gronwall's lemma. On the other hand, in the particular case where 
the initial vorticity is a single Dirac mass (of arbitrary size), uniqueness of the solution of (jl.2p 
was proved recently by C.E. Wayne and the author [19], using a dynamical system approach. 
An alternative proof of the same result can also be found in [16]. Finally, in the general case, 
it is possible to obtain uniqueness of the solution of ()1.2p by isolating the large Dirac masses in 
the initial measure and combining the approaches of [21] and |19j . This last step in the proof of 
Theorem 11.11 was achieved by I. Gallagher and the author in [15j . 

When the initial vorticity n = ad is a multiple of the Dirac mass (located at the origin) , the 
unique solution of p.2p is an explicit self-similar solution called the Lamb-Oseen vortex: 

a „/ X \ , G 



where 

The circulation parameter a G R measures the intensity of the vortex. It coincides, in this 
particular case, with the integral of the vorticity oj over the whole plane R^, a quantity which 
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is preserved under evolution. The dimensionless quantity \oi\lv is usually called the circulation 
Reynolds number. 

According to Theorem ll.il the vorticity equation (jl.2p has a unique global solution for any 
initial measure and any value of the viscosity parameter. It is very natural to investigate the 
behavior of that solution in the vanishing viscosity limit, especially if the initial data contain 
non-smooth structures such as point vortices, vortex sheets, or vortex patches. Indeed, if the 
viscosity is small, these structures will persist over a sufficiently long time scale to be observed 
and to influence the dynamics of the system. This question, however, is very difficult in its full 
generality, because Eq. ()1.2p reduces formally, as i/ — > 0, to the Eulerian vorticity equation 



which is not known to be well-posed in such a large space as M{M.'^). If the initial vorticity 
fj. S M.{M?) belongs to H^^{M?), and if the singular part of fi has a definite sign, then Eq. (jl.Sp 
has at least a global weak solution [131135] . but this result does not cover the case of point vortices 
due to the assumption fi £ H~^(M?). Furthermore, if we want to prove that the solution of (jl.Sp 
is unique, we have to assume that the initial vorticity is bounded [52] or almost bounded [53|[5T]. 

From a more general point of view, it is relatively easy to show that solutions of the Navier- 
Stokes equations converge, in the vanishing viscosity limit, to solutions of the Euler equations 
if we restrict ourselves to smooth solutions in a domain without boundary [14^ H8| [26| [2] . The 
situation is completely different in the presence of boundaries, especially if one chooses the 
classical non-slip boundary conditions for the solutions of the Navier-Stokes equations. In that 
case Prandtl boundary layers may form, and the inviscid limit becomes an extremely difficult 
question which has been rigorously treated so far only for well-prepared analytic data [46| [22]. 
But even in the absence of boundaries, hard problems can occur in the inviscid limit if one 
considers non-smooth solutions. Classical examples are listed below, in increasing order of 
singularity: 

1) Vortex patches. The simplest example in this category is the case where the vorticity is the 
characteristic function of a smooth bounded domain in R^. The corresponding velocity field is 
Lipschitz continuous, or almost Lipschitz if the boundary of the patch has singularities. The 
first convergence results, due to P. Constantin and J. Wu [318], and to J.-Y. Chemin [5j, hold 
for a general class of solutions including two-dimensional vortex patches. Several improvements 
have been subsequently obtained, especially by R. Danchin [11[ 112], H. Abidi and R. Danchin 
m, T. Hmidi [24[ 125]. and N. Masmoudi [41]. Closer to the spirit of the present work, we also 
quote a recent paper by F. Sueur [47j, where viscous transition profiles at the boundary of a 
vortex patch are systematically constructed, and provide a complete asymptotic expansion of 
the solution in powers of (z^t)-^/^. 

2) Vortex sheets. In a two-dimensional setting, this is the case where the initial vorticity is 
concentrated on a piece of curve in M^. The tangential component of the velocity field is 
discontinuous along the curve, thus creating a shear flow which is responsible for the celebrated 
Kelvin-Helmholtz instability. Due to this underlying instability, the inviscid limit for vortex 
sheets is at least as difficult as for Prandtl boundary layers, and has been rigorously treated so 
far only in the case of well-prepared analytic data [BJ [3] . 

3) Point vortices. This is the most singular example in our list, since here the velocity field is 
not even bounded near the vortex centers. However, in the case of point vortices, the evolution 
of the inviscid solution is given, at least formally, by the Helmholtz-Kirchhoff system which 
does not exhibit any dynamical instability. Therefore, one can reasonably hope to control the 
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vanishing viscosity limit in this particular situation. The first rigorous results in this direction 
were obtained by C. Marchioro [371 138j . and the aim of the present paper is to show how these 
results can be complemented to obtain an accurate description of the slightly viscous solution 
of (|1.2p when the initial condition is a finite collection of point vortices. 

After this general introduction, we now give a precise definition of the problem we want to 
study, and we state our main results. 

2 The viscous A'^-vortex solution 

Let be a positive integer. We take xi, . . . , xj\f G M? such that Xi 7^ Xj for i 7^ j, and we also 
fix ai, . . . , On G M \ {0}. Given any 1/ > 0, we denote by u;'^{x,t), u^{x,t) the unique solution 
of the vorticity equation (jl.2p with initial data 

N 

fi = ^ai(5(- - Xi) . (2.1) 

1=1 

This initial measure, which does not depend on the viscosity v, describes a superposition of N 
point vortices of circulations ai, . . . , located at the points xi, . . . , xn in R^. Our goal is to 
describe the behavior of the vorticity uj'^{x,t) in the vanishing viscosity limit. 

The measure fi is very singular, and we do not know how to construct even a weak solution 
of Euler's equation (II. 8p with such initial data. However, the inviscid motion of point vortices in 
the plane has been investigated by many authors, starting with H. von Helmholtz [23j and G. R. 
KirchhofF [28] who derived a system of ordinary differential equations describing the motion of 
the vortex centers. It is therefore reasonable to expect, in our case, that 

N 

uj''{-,t) ^ J^'^i^i- - ^ asz^^O, (2.2) 

i=l 

where z{t) = (zi{t), . . . , Z]\f{t)) denotes the solution of the Helmholtz-Kirchhoff system 

-^(*) = ^ \z,it) - z,it)\' ' "'^ ^ ^ • ^ ^ ^ 

As was already mentioned, the expression in the right-hand side of (j2.2p is not a weak solution 
of Euler's equation, because in that case the self-interaction terms in the nonlinearity u ■ Vuj 
are too singular to make sense even as distributions. However, as was shown by C. Marchioro 
and M. Pulvirenti \36\ [39], it is possible to derive system ()2.3p from Euler's equation by a 
rigorous procedure, which consists in approximating the Dirac masses in the initial data by 
small vortex patches, of size e > and circulations ai, . . . , a^, whose centers are located at 
the points xi, . . . , xat. Then the corresponding solution of ()1.8p converges weakly, as e — > 0, to 
the expression ()2.2p where the vortex positions zi{t), . . . , zj\f{t) are solutions of p.3p . We also 
recall that system (j2.3p is not globally well-posed for all initial data, because if A > 3 and if 
the circulations a, are not all the same sign, vortex collisions can occur in finite time for some 
exceptional initial configurations [40] . 

Our first result shows that the solution uj'^(x,t) of Eq. (II. 2p given by Theorem [LTj con- 
verges weakly, in the vanishing viscosity limit, to a superposition of point vortices which evolve 
according to ()2.3p . provided that vortex collisions do not occur. 
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Theorem 2.1 Assume that the point vortex system (j2.3p is well-posed on the time interval 
[0,T]. Then the solution uj'^{x,t) of the Navier-Sokes equation (II. 2p with initial data (12. ID 
satisfies 



N 

uj''{;t) . Va,5(--z,(t)) , forallte[0,T], (2.4) 

^ — ' 

i=l 

where z{t) = {zi{t), . . . , ZAr(t)) is the solution of (j2.3p . 

This theorem is closely related to a result by C. Marchioro [37^ I38j . which we now briefly 
describe. Instead of point vortices, Marchioro considers initial data of the form 

N 



i=l 

where, for each i € {1, . . . , A'^}, cof is a smooth vortex patch with a definite sign, which is centered 
at point Xi G M^, has compact support of size 0{e), and satisfies 

uji{x)dx = a1 > Qj . 

Under these assumptions, it is proved that the solution uj'^''^{x,t) of (|1.2p with initial data uJq 
converges to the expression in the right-hand side of (j2.2p in the double limit — > 0, e — > 0, 
provided 

< t'o 5 for some I'o > and /3 > . (2.5) 

Theorem 12.11 above corresponds to the limiting case e = 0, — > 0, which is precisely excluded 
by hypothesis (j2.5p . It should be mentioned, however, that restriction (|2.5p can be removed if 
the circulations ai, . . . ,aN are all the same sign, in which case Theorem 12.11 mav probably be 
established using the techniques developped in [37] . 

Marchioro's proof is based on a decomposition of the solution a;'^'^(x,f) into a sum of 
viscous vortex patches. The main idea is to control the spread of each patch by computing its 
moment of inertia with respect to a suitable point z^'^it), which is an approximate solution of 
(j2.3p . This argument does not give any information on the actual shape of the vortex patches, 
and is therefore not sufficient to provide a qualitative description of the solution uj^''^{x,t) for 
small e, u. Theorem 12.11 above suffers exactly from the same drawback. Its main interest is to 
provide a natural and rigorous derivation of the point vortex system ()2.3p . which differs from 
the classical approach of [36l [39] . 

The main goal of the present paper is to obtain a quantitative version of Theorem 12.11 which 
specifies the convergence rate in (12. 4p and provides a precise asymptotic expansion of the A^- 
vortex solution lij'^(x, t) in the vanishing viscosity limit. As in Marchioro's approach, our starting 
point is a decomposition of uj'^{x,t) into a sum of viscous vortex patches, which correspond 
to the atoms of the initial measure (12.11) . 



Lemma 2.2 The N-vortex solution can he decomposed as 

N N 

io''{x,t) = Y.ujUx,t) , ^x^(x,t) = ^uUx,t) , (2.6) 

1=1 i=l 

where, for each i £ {1,...,A^}, G C°((0, oo), ^^(M^) nL°°{R'^)) is the (unique) solution of 
the convection- diffusion equation 

-^ + (u''-VM = i.Aoj!^ , with co-^{.,t) ^ a,6i--x,) , (2.7) 
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and the velocity field is obtained from u;^ via the Biot-Savart law (jl.3p . Moreover LO^{x,t) 
has the same sign as ai for all x G and all t > 0, and f^2^i{x,t) dx = ai for all i > 0. 
Finally, there exists Kq > (depending only on v and |a| = |ai|+--- + |aAr|j such that 

K(..«)l<A-„^exp(-^), (2.8) 

for all i £ {1, ... , N}, all x G M^, and all t > 0. 

The proof of Lemma 12.21 which is borrowed from [19\ I15j . wih be reproduced in Section [5] 
for the reader's convenience. For the time being, we just point out that estimate (j2.8p gives a 
precise information on the A^-vortex solution in the hmit t — > for any fixed u > 0, but cannot 
be used to control in the limit ^ for fixed t > 0, because the constant Kq{i', \a\) 

blows up rapidly as u ^ 0. 

In the very particular case where = 1, we know from [19tll6j that ij'^{x, t) is just a suitable 
translate of Oseen's vortex (jl.6p . From now on, we assume that N > 2, and we suppose that 
the point vortex system ()2.3p is well-posed on the time interval [0, T]. We denote 

d = min min \zi(t) - Zj(t)\ > , (2.9) 
ie[0,T] i^i ' ' ' " 

and we also introduce the turnover time 

To = - — - , where lal = |ai| + • • • + \ai\[\ ■ (2-10) 
|a| 

Our goal is to show that, if > is sufficiently small, the A^-vortex solution uj'^{x,t) looks like 
a superposition of Oseen vortices located at some points Zi(t), . . . ,zff(t) which satisfy the 
following viscous regularization of the Helmholtz-Kirchhoff system: 

,"(0)^.,. (2.n) 



dt ^ ^/JJi V ^ 

where v'^ is given by (jl.7p . The reason for using that system instead of (j2.3p will be explained 
at the beginning of Section [3l For the moment, we just observe that system (12. lip is globally 
well-posed for positive times, and that the solutions (t) are exponentially close to the solutions 
Zi{t) of ()2.3p if the viscosity is sufficiently small. 

Lemma 2.3 Assuming pairwise distinct initial positions xi, . . . ,xn, system (j2.1ip is globally 
well-posed for positive times, for any value of u > 0. Moreover, if the solution of (j2.3p satisfies 
(j2.9p . there exists Ki > (depending on the ratio T/Tq) such that 

1 / \ 

- max \z'i{t) - Zi{t)\ < Kiey.Y>[-—] , for all t e {0,T] . (2.12) 

d i=l,...,N \ OVt/ 

The proof of Lemma 12.31 is also postponed to Section [5l We refer to [Uj for a discussion of the 
validity of system (j2.1ip as a model for the dynamics of interacting vortices. 

To obtain a precise description of each vortex patch uj'^{x,t) in a neighborhood of z^{t), we 
introduce, for each z G {1, . . . , A^}, the self-similar variable 

^ _ x-zUt) 



ut 
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and we define rescaled functions 'w'^{S,,t) £ M and t) E M? by setting 



Given a small /? G (0, 1), which will be specified later, we also introduce a weighted space X 
defined by the following norm: 



w{0\' e^^^^/U^y^^ . (2.14) 



We already know from [15] that tff t) converges to the Gaussian profile G(^) as t ^ for any 
fixed > 0. Our first main result shows that a similar result holds in the vanishing viscosity 
limit, uniformly in time on the interval (0, T]. 

Theorem 2.4 Assume that the point vortex system (j2.3p is well-posed on the time interval 
[0,T], and let oj"{x^t) he the solution of (|1.2p with initial data (j2.1|) . Ifuj'^{x,t) is decomposed 
as in ()2.6p . then the rescaled profiles t) defined by (j2.13p satisfy 

max Ww'ii-^t) -G\\x = o(%^ , as z/ ^ , (2.15) 



uniformly for t G (0,T], where d is given by (I2.9p . 

More precisely, the proof shows that there exist positive constants P and K2, depending on 
the ratio T/Tq, such that 

max |K(-,t)-G||x < -f^2 ^ , (2.16) 

i=l,...,N 

for all t £ (0, T], provided u is small enoug h so that K2{vT/d'^) < 1. This result means that, 
when the viscosity u is small, the A'^-vortex solution uj^{x^t) looks like a superposition of N 
Oseen vortices located at the points Zj'(t), . . . , zf^{t), which evolve in time according to (|2.1ip . 
Since X ^ L^(M^) and since z'j^{t) Zi{t) as 1/ —> by Lemma [231 estimate ()2.15p implies in 
particular (j2.4p , and Theorem 12.11 is thus a direct consequence of Theorem 12.41 

Theorem 12.41 is already very satisfactory, but it does not seem possible to prove it directly 
without computing a higher order approximation of the A^- vortex solution. This rather surprising 
claim will be justified in Section [3l but for the moment it can be roughly explained as follows. 
As is clear from ()1.6p . the velocity field of Oseen's vortex is very large near the center if the 
viscosity v is small, with a maximal angular speed of the order of |a|/(z/t). As long as the 
vortex stays isolated, it does not feel at all the effect of its own velocity field, because of radial 
symmetry. However, if a vortex is advected by a non-homogenous external field, which in our 
case is produced by the other A'^— 1 vortices, it will get deformed and, consequently, will start 
feeling the influence of its own velocity field. If the ratio {al/u is large, this self-interaction will 
have a very strong effect, even if the deformation is quite small. In particular, one may fear 
that the vortex gets further deformed, increasing in turn the self-interaction itself, and that 
the whole process results in a violent instability. In fact, this catastrophic scenario does not 
happen. Remarkably enough, a rapidly rotating Oseen vortex in an external field adapts its 
shape in such a way that the self-interaction counterbalances the strain of the external field 
[50\ [59] . This fundamental observation will be the basis for our analysis in Section [3l It explains 
why one can observe, in turbulent two-dimensional flows, stable asymmetric vortices which in 
a first approximation are simply advected by the main stream. The same mechanism accounts 
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for the existence of stable asymmetric Burgers vortices, which are stationary solutions of the 
three-dimensional Navier-Stokes equations in a linear strain field [431 \20\ [33| . 

To compute the self-interactions of the vortices, which play a crucial role in the convergence 
proof, we construct a higher order approximation of the A'^-vortex solution in the following way. 
For each i £ {1, . . . , N} and all t G [0, T], we denote 

y^ntt) = G{o + Q){m,t)+F,^{tt)} , C€R\ (2.17) 

where Fi{^,t) is a radially symmetric function of whose precise expression is given in (I3.42p 
below, and F^[^^t) is a nonsymmetric correction which satisfies 



where Zij{t) = z^(t) — Zj{t). Here w : (0, oo) — > M is a smooth, positive function satisfying 

uj{r) Cir^ as r ^ and uj{r) (72r^e~'"^/^ as r — > oo for some Ci,C2 > 0, see Eq. ()3.28p 
below. The right-hand side of (12.17P is the beginning of an asymptotic expansion of the rescaled 
vortex patch w^{S,,t) in powers of the non-dimensional parameter {iyt)/d?. Each term in this 
expansion can in turn be developped in powers of z^/|a|. The most important physical effect is 
due to the nonsymmetric term F!^{S^, t), which describes to leading order the deformation of the 
i*^ vortex due to the influence of the other vortices. Keeping only that term and using polar 
coordinates ^ = {r cos 9, r sin 9), we can rewrite (|2.17p in the following simplified form 

= 9ir) + '^Y.^^-^cos{2{9-9,,{t))) + ... , (2.19) 

where g{\(,\) = G(^) and %(t) is the argument of the planar vector Zij(t) = z^(t) — Zj{t). This 
formula allows to compute the principal axes and the eccentricities of the vorticity contours, 
which are elliptical at this level of approximation. 

Using these notations, our final result can now be stated as follows: 

Theorem 2.5 Assume that the point vortex system (j2.3p is well-posed on the time interval 
[0,T], and letuj'^{x,t) be the solution of ()1.2p with initial data ()2.ip . Ifuj'^{x,t) is decomposed 
as in (|2.6p . then the rescaled profiles w'^{S,,t) defined by (j2.13p satisfy 

.max^||<(.,t)-u;r(-,t)|U = , as u ^ , (2.20) 

uniformly for t £ (0,T], where w^"^"^ is given by (j2.17p . 

As is clear from (|2.15p . (j2.17p . (j2.20p . Theorem 12.51 implies immediately Theorem 12.41 hence 
also Theorem 12.11 Note that the convergence result is now accurate enough so that the first 
order corrections to the Gaussian profile in (12.17P are much larger, for small viscosities, than 
the remainder terms which are summarized in the right-hand side of (j2.20p . This means in 
particular that the deformations of the interacting vortices are really given, to leading order, by 
()2.19p . According to that formula, each vortex adapts its shape instantaneously to the relative 
positions of the other vortices, without oscillations or inertia. In this sense, point vortices can 
be considered as well-prepared initial data, and the A'^-vortex solution is an example of the 
"metastable regime" described in the introduction. In contrast, one should mention that the 



9 



first order radially symmetric corrections t) do not only depend on the instantaneous vortex 
positions, but on the whole history of the system, see (I3.42p . 

The rest of this paper is devoted to the proof of Theorem 12.51 which is divided in two main 
steps. In Section [3l we construct an approximate solution of our system, with the property 
that the associated residuum is extremely small in the vanishing viscosity limit. This approxi- 
mation differs from (j2.17p by higher order corrections which are necessary to reach the desired 
accuracy, but will eventually be absorbed in the right-hand side of (j2.20p . As is explained in 
Section 13.31 a difficulty in this construction comes from the fact that the radially symmetric 
and the nonsymmetric terms in the approximate solution w^^^{^,t) have a different origin, and 
play a different role. Once a suitable approximation has been constructed, our final task is to 
control the remainder ■Wi{S^,t) — w^'^^{(,,t) uniformly in time and in the viscosity parameter u. 
This will be done in Section HI using appropriate energy estimates. A major technical difficulty 
comes here from the fact that we do not want to assume that the ratio T/Tq is small. If we did 
so, the proof would be considerably simpler, and we could replace the weig 
in the definition (j2.14p of our function space X, thus improving our convergence result. In the 
general case, however, we have to use a rather delicate energy estimate involving time-dependent 
weights Pi{S,,t), which will be constructed in Section [4.11 The price to pay is a slightly weaker 
control of the remainder, and the fact that all our constants, such as K2 and (3 in (j2.16p . have 
a bad dependence on T if T Tq. 

In conclusion, our results show that the vanishing viscosity limit can be rigorously controlled 
in the particular case of point vortices, due to the remarkable dynamic and strucural stability 
properties of the Oseen vortices. These properties, which were established in [El |20l |32] , play a 
crucial role both in the construction of the approximate solution in Section O and in the energy 
estimates of Section HI 



3 Construction of an approximate solution 

In this section, we show how to construct an asymptotic expansion of the A^-vortex solution in 
the vanishing viscosity limit. Our starting point is the evolution system satisfied by the rescaled 
profiles w^{S,,t), which from now on will be denoted by ■Wi{^,t), Vi{^,t) for simplicity. 

Replacing (|2.13p into (j2.7p . we obtain, for i = 1, . . . , A^, 

t5t«;.(C,t) + |^^^t;,(c + ^,t)-y^zKi)|-V«;.(e,t) = iCw.mt), (3.1) 

where 

Cw = Aw + -C-Vw + w . (3.2) 

Here and in what follows, we denote the vortex positions by Zi{t) instead of (t), to keep the 
formulas as simple as possible. We also recall that Zij{t) = Zi{t) — Zj{t). 

The initial value problem for system (j3.ip at time t = is not well-posed, because the 
time derivative appears in the singular form tdf. A convenient way to avoid this difficulty is to 
introduce a new variable r = log(t/T), so that dr = tdt- With this parametrization, the solution 
of (jS.ip given by Lemma 12.21 is defined for all r G (— oo,0], and converges to the profile G of 
Oseen's vortex as r ^ —00, see |15[ Proposition 4.5]. For simplicity, we keep here the original 
time t, because this is the natural variable for the ODE system (j2.3p or (j2.1ip . 
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3.1 Residuum of the naive approximation 



If w{^,t) = {wi{(,,t), . . . ,WN{C,t)) is an approximate solution of system (13. ip . we define the 
residuum t) = t), . . . , Rn{S,j t)) of this approximation by 

N 

n ■ / r • • I r 1 \ It , 



^ 1=1 ^ 



for alH G {1, . . . , A^}. Here and in the sequel, it is always understood that Vi{^, t) is the velocity 
field corresponding to Wi{^,t) via the Biot-Savart law ()1.3p . 

In view of Theorem l2.4l the solution of (jS.ip we are interested in satisfies Wi{£,, t) ~ G(^) and 
'ViiC^i) ~ ^'^(0 for all t G (0, T], if u is sufficiently small. Since dtG = CG = 0, the residuum of 
this naive approximation (where Wi(^, t) = G(^) for alH € {1, . . . , A^}) is 

This expression looks singular in the limit ^ 0, but the problem can be eliminated by an 
appropriate choice of the vortex positions 2:i(t), . . . , ziy{t). Indeed, in (jS.Sp . the quantity inside 
the curly brackets {•} vanishes for ,^ = if we set 



N 



m-Y^^i^), (3.4) 

This is exactly the regularized point vortex system, which was already introduced in (j2.1ip . 

Prom now on we always assume, as in Theorems 12.41 and 12. 5( that the original point vortex 
system ()2.3p is well-posed on the time interval [0, T]. For any > 0, we denote by 2:i(t), . . . , ZN{t) 
the solution of ()3.4p with initial data xi, . . . ,X7v- In view of Lemma 12.3^ we can assume that 
this solution satisfies (j2.9p for some d > (independent of v), provided v is sufficiently small. 
Replacing (j3.4p into (|3.3p , we obtain the following expression of the residuum 



«f«..)^E?{"°(«-^)-.'1^)}-VO(0. (3.) 

Remark that the sum now runs on the indices j ^ i, because the term corresponding to j = i 
vanishes. Our first task is to compute an asymptotic expansion of the right-hand side of ()3.5p 
in the vanishing viscosity limit. 

Proposition 3.1 For i = 1, . . . , N , we have 



Rri^^t) = -^[MC^t) + [^) i?,(^,t) + (^-jQ(e,t) + i?r(C,i)| , (3.6) 
for all ^ G and all t £ (0, T], where 

27r ^ Fii(*) 
^^(^'*) = £ E g ^t'f r - 4(e • -m') GiO , (3.7) 

«(f . ') = ? E f '^■-•^''>'«;;-''"'' (2« . - itf GK) . 
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Moreover, for any 7 < 1, there exists C > such that 



l/t\3/2 



-7l«lV4 



< t < T . 



(3.8) 



Remarks. Proposition 13.11 provides an expansion of the residuum F6^\^,t) in powers of the 
dimensionless parameter (i/t)^/^/d, where (z^t)^/^ is a diffusion length which gives the typical size 
of each vortex at time t, and d is the minimal distance between the vortex centers, see ()2.9p . As 
is clear from (j3.7p . the parameter d in (|3.6p has been introduced rather artificially, to ensure that 
the quantities j4j, Bi, Ci are dimensionless, as is the residuum itself. The proof will show that 
an expansion of the form (j3.6p can be performed to arbitrarily high orders, but for simplicity we 
keep only the terms which will be necessary to prove Theorems 12.41 and 12.51 Finally, we remark 
that the prefactor Oit/d"^ in (j3.6p is bounded by t/To, where Tq is the turnover time introduced 
in (pro]) . 

Proof. Fix 7 G (0, 1), and let 71 = 1 — 7. Since VG(^) = — |^G(^) decreases rapidly as |^| — > 00, 
it is clear that the residuum (|3.5p is extremely small if |^| is large. For instance, if |^| > d/ {2\/id), 
we can bound 



\R^\c,t)\ <--iiv 



a 



«"-^ieiG'(o 



< c 



a 



t /(i2\3/2 



d^ \utJ 



- exp 



1 |a| 
47r V 



I6ut 



) 



-7I5IV4 



-7i|?lV4 g-7l«lV4 



where in the last inequality we have used the fact that \^\e "^^l^l^/^ is a decreasing function of 
|,^| when 1^1 ^ 1. A similar argument shows that 

\M^,t)\ + mc,t)\ + \c\iu)\ < ^(S'''''p("i^)"~^'^''^' ' 

if Id ^ dl(2\pui). We conclude that expansion (|3.6p holds in that region, with a remainder term 
satisfying a much better estimate than (|3.8p . 

We now consider the case where 1^ < d/{2\/id). Since |^;jj(t)| = \zi{t) — Zj{t)\ > dhy (j2.9p 
if z 7^ we have 



Zij{t) 



> 



d 



> 2|d , and 



> 



d 



2Jvt 



(3.9) 



To estimate the right-hand side of (j3.5p . we have to compute the difference v'^{^ + v) ~ '^'^{'n)-, 
with r] = Zij{t) / ^/vi. Using definition (jl.7p . we obtain the identity 



v^{i + rii)-v'^{vi) 



where 



2tt 



(3.10) 



rr 



-|r,|V4 



-|C+r,|2/4 



In particular, it follows from (j3.9p that 



Zij{t) 



exp 



d' \ 



IQut 



) 



(3.11) 



hence the contributions of the term V2 to the residuum (|3.5p are exponentially small and can be 
incorporated in the remainder term. To compute we use the following elementary lemma: 
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Lemma 3.2 For all ^,77 G with \^\ < \ri\, we have 



°° I tin 

?-^i(e,r/) = 5](-ir-i^sm(n(0-0)) , (3.12) 

n=2 

where 9 denotes the polar argument of ^ and (p the argument of ij. 

For completeness, the proof of Lemma 13.21 will be given in Section [5l Applying (|3.12p with 
r/ = Zij{t) / ^/vi and keeping only the first three terms in the expansion, we obtain 

J . («. ^) = - ,m(2(« - *)) + raW! _ ^, 

+ (3.13, 

I ^ij I ^ I ^ij I ^ 

where — </) is the signed angle between ^ and Zij. In particular, we have the relations 



sin(0 -(/>) = , cos(0 - 4>) 



\^\\Zij\ ' \^\\Zij\ ' 

in terms of which the higher order trigonometric expressions appearing in ()3.13p can be computed 
using the well-known formulas sin(2a) = 2 sin(a) cos(a), sin(3a) = sin(a)(4cos2(a) — 1), and 
sin(4a) = 4sin(a) cos(a)(2cos2(a) — 1). Summarizing, we have shown that 



- - / 4vr ^ . V V-t 

= ^[Aitt) + QY^'m^t) + {^)c^{^,t) + MC,t) 

for |,^| < d/{2y/iA), where ^j, Bi, Ci are given by (13. 7p and Ri{S,,t) satisfies the bound (j3.8p . As 
was already observed, the same result holds if we replace Vi with Vi + V2, so we conclude that 
expansion (13. 6p is valid in the region |^| < d/{2\/lA) too. The proof of Proposition 13.11 is thus 
complete. □ 

It is important to remark that the residuum Rf^\^,t) does not converge to zero as — > 0, 
because of the leading order term Ai{^,t). If we decompose the solution of (13. ip as Wi{S,,t) = 
G(^) + Wi{S,j t)j equation for Wi{^, t) will contain a source term of size 0{1) as ^ 0, and we 
therefore expect that the remainder Wi{(,,t) itself will be of size 0{1) after a short time. But, 
as is easily verified, the equation for t<;j(^,t) contains nonlinear terms with a prefactor of size 
0(1/"^), and such terms cannot be controlled in the vanishing viscosity limit if t(;j(^,t) is 0{1). 
This is the reason why it is necessary to construct a more precise approximate solution of (13. ip . 
with a sufficiently small residuum, in order to desingularize the equation for the remainder. 



3.2 First order approximation 

We look for an approximate solution of (13. ip of the form 

where, for each i G {1,...,A^}, Fi{S^,t) is a smooth vorticity profile to be determined, and 
i;^'(^,i) is the velocity field obtained from Fi(^,t) via the Biot-Savart law (|1.3p . In fact, we 
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shall need later a more precise approximation of the A^-vortex solution, but we prefer starting 
with ()3.14p to describe the procedure in a relatively simple setting. Our goal is to chose the 
profile -Fj(^,t) so as to minimize the residuum of our approximation, which is 



N 



V 



Zij{t) 



Ut 



Zij{t) 



Using (IXT4D and the definition of i?f 

we find after some calculations 



?{0)/ 



ad 



where 



(3.15) 



VFi{tt)+v''^{^ + -^^,t)-VG{0 



... 

It is easy to check, at least formally, that R^^\tt) is Oiut/d?) as — > (this calculation will 
be done rigorously later, when the profile Fi will be determined). So it follows from (|3.6|) and 
(f3l^ that 



(3.16) 



To minimize the residuum, it is natural to impose Ai + v^- VFi + v^' • VG = 0. As we shall see, 
this "elliptic" equation has a solution, but does not completely determine the profile Fi. 
To prove this claim, we first introduce some notation. Let Y denote the Hilbert space 



equipped with the scalar product 

{wi,W2)y 



2m2\ 



/ mc)pei«i'/4de < oo > , 

Jr2 J 



(e)«^2(e)ei«i'/^de . 



(3.17) 



(3.18) 



We consider the linear operator A : D{A) Y defined by -D(A) = {w £ Y \ v'~^ ■ Vw G 1"} and 



Aw = - Vw + v -VG 



w G D{A) 



(3.19) 



where (as always) v denotes the velocity field obtained from w via the Biot-Savart law (jl.3p . 
With these notations, the equation we have to solve becomes AFi + Ai = 0, and we are therefore 
interested in computing the (partial) inverse of A on appropriate subspaces. This operator has 
been extensively studied, because it plays a prominent role in the stability properties of the 
Oseen vortices and the construction of asymmetric Burgers vortices [191 l20| [32| [33\ I34j . In 
particular, we have 



-A. Moreover, 
(3.20) 



Proposition 3.3 /7g|. [3^ The operator A is skew-adjoint in Y , so that A* 

Ker(A) = Yq® {/Ji^iG + /^s^sG | /32 G M} , 
where Yq <ZY is the subspace of all radially symmetric functions. 
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If w £ Yq, the corresponding velocity field v satisfies C ' = 0, hence (j3.19p immediately 
implies that Aw = 0. On the other hand, if we differentiate the identity v'^ ■ VG = with respect 
to ^1 and ^2, we obtain A{diG) = A{d2G) = 0. So the right-hand side of ()3.20p is certainly 
contained in the kernel of A, and the converse inclusion was proved in [32]. On the other hand, 
the fact that A is skew-adjoint in Y implies that Ker(A) = Ran(A)-'-, so we know that the range 
of A is a dense subspace of Ker(A)"'". We shall now prove that Ran(A) contains Ker(A)-'" n Z, 
where 

Z = : ^ M el^l'/^u- e ^(M^)} C Y . 

Here 5(M^) denotes the space of all smooth, rapidly decreasing functions on M?. 

As was observed e.g. in [19], the operator A commutes with the group 5*0(2) of all rotations 
about the origin. It is thus natural to decompose 



oo oo 
Y = Yn = PnY 

n=0 n=0 



where P„ is the orthogonal projection defined in polar coordinates (r, 6) by the formula 

2 J /"Ztt 

{Pnw){r COS 6 ,r sulO) = ^ / w{r cos 6' ,r sin 6') cos{n{6 — 6')) dO' , nGN. 

Then Yq = PqY is the subspace of all radially symmetric functions, and for n > 1 the subspace 
Yn = PnY contains all functions of the form w{r cos 9 ,r sm. 9) = ai{r)cos{n9) + a2ir) s'm{n9). 
It is not difficult to verify that the projections P„ commute with A for all n G N, and explicit 
formulas for the restrictions A„ = P„AP„ are given in \19\ Section 4.1.1]. To formulate the main 
technical result of this section, we shall use the following notations: 

^„ = _±_(i_,-.v4,, ,(,) = |!i . _!;!Zi_ , ,> 0.(3,21) 

Lemma 3.4 If z GYnCiZ for some n>2, there exists a unique w GYnCi Z such that Aw = z. 
In particular, if z = a{r) sm{n9), then w = —uj{r) cos{n9), where 

uj{r) = h{r)Q{r) + , r > , (3.22) 

n(p[r) 

and $7 : (0, oo) is the unique solution of the differential equation 

_ i(rl7'(r))' + (4 - h{r))nir) = ^ , r > , (3.23) 
such that Q.[r) = 0{r^) as r ^ and il(r) = ©(r^") as r ^ oo. 



Remarks. 

1. By rotation invariance, if z = a{r)cos{n9), then w = uj{r) sin{n9), and the relation between 
uj and a is unchanged. The general case where z = ai{r)cos{n9) + a2{r)sin{n9) follows by 
linearity. 

2. The conclusion of Lemma 13.41 is wrong for n = 1. Indeed, since dkG = -^CkG for A; = 1,2, 
it is clear that diG,d2G £ YiCi Z, but Proposition 13.31 asserts that these functions belong to 
Ker(A) = Ran(A)-'". However, if z G Yi n Z satisfies {z,dkG)Y = for k = 1,2, one can show 
that there exists a unique w G Yi f] Z n Ker(A)-'- such that Aw = z. This result will not be 
needed in what follows, so we omit the proof. 
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3. If w is as in Lemma 13.41 and if v is the velocity field associated to w via the Biot-Savart law 
(ll.3p . the proof will show that v is smooth and satisfies 



\vm = 0{\(r') as e - , and \viO\ = 0(|Cr"~') as |e| - oo . (3.24) 

Proof of Lemma 13.41 A particular case of Lemma 13.41 was proved in [20' Proposition 3.1]. 
Since the general case is quite similar, we just indicate here the main steps and refer to [2D] for 
further details. 



Assume that z £ YnOZ for some n >2. By Proposition [3i3l we have z £ Ker(A) = Ran(A). 
Our task is to verify that z G Ran(A), and that there exists a unique w £ Yn D Z such that 
Aw = z. Without loss of generality, we assume that z = a{r)sm{n6). Then a : ^ M is a 
smooth function with the property that /^a(r) decays rapidly as r — > oo. Furthermore, we 
can write a(r) = r"yl(r^), where A : [0, oo) — > R is again a smooth function. In particular, 
we have a(r) = 0{r'^) as r — > 0. We look for a solution w of the form w = —uj{r)cos{n6). 
The corresponding stream function, which is defined by the relation — A^f = w, satisfies ^ = 
—Q{r)cos{n9), where is the unique regular solution of the differential equation 

- -{rQ'ir))' + -2 n{r) = co{r) , r > . (3.25) 
Moreover, the velocity field v = — V"'"^' has the following expression 

71 

V = — ^}{r) sin{n6) er + ^l' (r) cos{n9) 60 , 
r 

where e,., eg are the unit vectors in the radial and azimuthal directions, respectively. Thus, 
using definitions (|1.7p and (j3.2ip . we obtain 

v-VG = ——Q.{r)g{r)s\\i{n9) , and ■'Vw = nuj{r)(l){r) sm{n9) . 

In particular, we see that Aw = z if and only if —^flg + rKpto = a, which is ()3.22p . Furthermore, 
combining (j3.22p and (j3.25p . we obtain the differential equation (j3.23p which determines O. 

It remains to verify that (|3.23p has indeed a unique solution with the desired properties. We 
first consider the homogeneous equation obtained by setting a(r) = in (j3.23p . Since /i(r) 1 
as r ^ and h{r) decays rapidly as r — > oo, this linear equation has two particular solutions 
ip+, V'- which satisfy 

'ip^{r) ~ r" as r — > , and ip+{r) ~ r"" as r ^ oo . (3.26) 

These solutions are of course unique. Moreover, since 77 > 2, the coefficient n^/r'^-h{r) in ([3:23|) 
is always positive, because sup^^g '^^^('") — 2.59... < 77^. It then follows from the Maximum 
Principle that the functions ip^, ip- are strictly positive and satisfy 

ip-{r) ~ K^r" as r — > oo , and ip+{r) ~ K+r"" as r — > , 

for some k_,k^ > 0. In particular, ip_ are linearly independent. In addition, their Wron- 
skian determinant satisfies W = — = wq/v for some wq > 0, and it follows that 

= K_ = wo/{2n). We now return to the full equation ()3.23p . and consider the particular 
solution given by the explicit formula 

n{r) = 7A+(r) rA^_(y)^dy + 7/;_(r) H ^^,^{y)^dy , r>0. (3.27) 
Jq wq ncply) Jr wo n(j){y) 
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As is easily verified, we have ^{r) = 0{r'^) as r ^ 0, ^{r) = 0{r~^) as r — > cxd, and Vl is the 
unique solution of ()3.23p with these properties, li w = —uj{r)cos{n9), where uj is defined by 
()3.22p . then Aw = z hy construction, and it is not difficult to see that w £ YnCi Z. Indeed, it is 
clear that w is smooth away from the origin, and the fact that e'^' ^^w decays rapidly at infinity 
follows immediately from (j3.22p . To prove that w is smooth in a neighborhood of zero, we 
observe that any regular solution of (I3.23P has the form 0(r) = r"<I>(r^), where : [0, oo) M 
is a smooth function. Using (j3.22p . we conclude that 

w{o = -\cr cosine) [hmm\') + ^^) 

is smooth also near the origin, because |^|"cos(n0) = Re((.^i + 1^2)^) is a homogeneous polyno- 
mial in ^. The proof of Lemma 13.41 is thus complete. □ 

Remark. Of course, the argument above fails if n = 1, because the coefficient — h{r) in 
(j3.23p is no longer positive. In fact, it is easy to verify that the functions ip+jip- defined by 
(|3.26p are linearly dependent in that case. 

Equipped with Lemma 13.41 we now go back to the determination of the vorticity profile 
Fi{^, t) in (|3.14p . The equation we have to solve is AFi + Ai = 0, where Ai{^, t) is given by ()3.7p . 
Using polar coordinates (r, 6) as before, we can write 

47r ^ ai \zij[t)\^ 

where Oij{t) is the argument of the vector Zij{t) = Zi{t) — Zj{t). This expression shows that 
Ai{-,t) G I2 n Z for any t € [0,T], hence by Lemma EH there exists a unique G I2 n Z 

such that AF° + = 0. Explicitly, 

F^{.,t) = ^^(r)5]^^-i-^ cos(2(0-0,,(t))) , (3.28) 

where Lo{r) is given by (j3.22p . (j3.23p with n = 2 and a(r) = r'^g{r). It follows in particular from 
(I222I), that 

\F?{i,t)\<c\i\\i + \i\^)e-\^\"i\ tG(o,r]. 

If we now return to (j3.15p and choose Fi{^,t) = F^{S,,t), it is easy to verify that 

so we succeded in constructing an approximate solution with a smaller residuum than uf'^ . As 
is explained in Section [21 the profile FP(^,t) describes to leading order the deformations of the 
vortices due to mutual interaction. 

Before going further, we state and prove a variant of Lemma 13.41 which will be useful in the 
next section. 

Lemma 3.5 Assume that z £ ¥„, f] Z for some n > 2, and let w € YnCi Z be the solution of 
Aw = z given by Lemma \3.4\ Then for all e ^ the equation 

€{1-C)w' + Aw^ = z (3.29) 
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has a unique solution ^YnD Z. Moreover, there exists C > (depending on z) such that 

\\w^ - w\\y < -^1^^ , foralle^Q. (3.30) 

Proof. Again, a particular case of Lemma 13.51 has been proved in |20t Proposition 3.4]. As is 
well-known (see e.g. |18]). the operator C defined by (j3.2p is self-adjoint in Y and its spectrum is 
given by (j{C) = {— § | re = 0, 1, 2, . . . }. Since A is skew-symmetric and relatively compact with 
respect to £, it follows that the operator —eC + A is maximal accretive for any e > 0. Thus, for 
any z ^Y, equation ()3.29p has a unique solution w'' G Y, which satisfies ||u)'^||y < e~"'^||z||y. A 
similar result holds of course for e < 0. 

We now consider the particular case where z G n Z for some re > 2. Since C commutes 
with the projection P^, it is clear that PnW^ = w'^, hence £ Yn- In that subspace, Eq. (I3.29P 
reduces to an ordinary differential equation which can be studied as in the proof of Lemma 13. 4[ 
In particular, it is straightforward to check that e'^' /^w'' decays rapidly as |^| oo, so that 

^Ynf] Z. On the other hand, since 

e(l - C){w' -w)+ A{w'' -w) = -e(l - C)w , (3.31) 

the fact that ||[e(l - £) + A]-i|| < |e|-i implies 

\\w^-w\\y < ||[e(l - £) + A]"^|| ||e(l - £)u'||y < ||(1 - £)u;||y , (3.32) 

hence \\w'' — w\\y is uniformly bounded for all e 7^ 0. 

Finally, since w & YnCi Z, we have (1 — C)w £ YnCi Z, hence by Lemma [33] there exists a 
unique w £ YnCi Z such that Aw = {1 — C)w. Then ()3.3ip takes the equivalent form 

[e(l -C) + A]{w^-w + ew) = e^{l - C)w , 

from which we deduce 

-w = -ew + e2[e(l - £) + A]"^(l - C)w . 
Since ||[e(l — C) + A]~^|| < |e|~-^, we conclude that 

\\w^ - w\\y < |e|(||w||y + 11(1 - C)w\\y) ■ (3.33) 
Combining (j3.32p and (|3.33p , we obtain (j3.30p . This concludes the proof of Lemma 13.51 □ 

Remark 3.6 For later use, we also introduce the space Z^ d Z defined by 




where 5=k(M^) denotes the space of all smooth functions u; : — > M such that w and all its 
derivatives have at most a polynomial growth at infinity. As is easily verified, Lemmas \3.4\ and 
13.51 still hold if we replace everywhere Z by Z^,. In particular, if z £ YnO Z^,, the solution of 
(j3.29p belongs to YnCi Z^. Thus, for any 7 < 1, there exists C > such that 

k'(OI + l^w'iOl < Ce-^l«l'/^ , ^£R\ 
and using for instance (I3.3ip one can show that the constant C is independent of e. 



18 



3.3 Third order approximation 

We now construct an approximate solution that will be accurate enough to prove Theorem 12. 51 
We set 

\a J \a J \a J ^^^^^^ 

where the vorticity profiles Fi, Hi, Ki have to be determined, and the velocity fields v^*, 
v^"- are obtained from Fi, Hi, Ki via the Biot-Savart law (11. 3p . As in Proposition 13.11 the main 
expansion parameter in (j3.34p is {vt/cP). The profiles Fi, Hi, Ki still depend on the viscosity 
i^, but they all have a finite limit as ^ 0. 

Our first task is to compute the residuum of the third-order expansion ()3.34p . as 

an approximate solution of (|3.ip . By a direct calculation, we find 



= [^){tdtFi + F,-£F?j{^,t) 



+ © (i^t + 1^^- CH^ {i, t) (3.35) 



+ (^) ' (tdtKi + 2Ki - CKi'^ it t) 

and using (j3.5p . (j3.19p we obtain 

ait f . ^ fi^i\h , -r-r f 



By definition, the residuum R- is the sum of all terms in (j3.35p and (|3.36p . 

Next, we separate the lower order terms, which will have to be eliminated by an appropriate 
choice of Fi,Hi,Ki, from the higher order terms, which will be automatically negligible. We 
thus decompose 

R?\^,t) = Rl{i,t)+R^{i,t) , 
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where collects the lower order terms in the residuum, namely 

Our goal is to choose the vorticity profiles Fj, Hi, Ki so as to minimize the quantity -Rf (C> t) in 
the vanishing viscosity limit. The contributions of the naive residuum ()3.6p are easily eliminated 
by successive applications of Lemma [3.41 we first take Fi{^, t) so that AFi + Ai = 0, then Hi{^, t) 
so that AHi + i?j = 0, and so on. In this way, we obtain a residuum Rf{^,t) of size 0{vt/cP), 
but unfortunately this is not sufficient to prove Theorems 12.41 and 12.51 To obtain a more precise 
estimate, we also need to eliminate all terms in the last two lines of ()3.37p . This requires a more 
careful choice of Fj, Hi, Ki, which will be done in three steps: 

1) First, we take 

F^{i,t) = m,t) + Fai,t) , e G , t e (o,r] , (3.38) 

where t) G Iq n Z and F^{-,t) E ^2 n Z. More precisely: 

a) The profile F^{^,t) is the unique solution of the elliptic equation 

— (Ff - CFn + AF,' + Ai = 0, (3.39) 

Ui 

as given by Lemma [331 Since Ai{-,t) G 1^2 n Z, we have F^{-,t) G I2 n Z for ah t G (0,r]. 
Moreover, by ([330]) . 

\\F^{;t)-Ff{.,t)\\y < , tG(0,r] , 

\OLi\ + V 

where F^[i,t) is given by (j3.28p . 

h) The profile Fi(^,t) is the unique solution of the linear parabolic equation 

tdtFi + Fi- CFi + ^{Po(^^^ • VF^) + PoiDi ■ VFn} = , (3.40) 

with initial data Fi{-,0) = 0. Here Pq is the orthogonal projection (in Y) onto the radially 
symmetric functions, and Di{^,t) is the divergence-free vector field given by 

It is clear that Fi{^,t) is well-defined. In fact, if S{t) = exp(r£) denotes the Co-semigroup in 
Y generated by £, we have the explicit formula 

,.f , 

H-.t) = S {log -y-P^Q,{;s)ds , (3.42) 

where Qi = • VF^ + A • Vi^'^. 
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2) Next, we determine the profile Hi{S,,t). From the proof of Proposition 13.11 it is clear that 
Bi{-,t) € Y^n Z. Thus, by Lemma 13.41 there exists a unique solution Hi{-,t) £ Y^f] Z to the 
equation 

AHii; t) + B,{; t) = 0, te{0,T] . (3.43) 

3) Finally, we set 

Ki{^, t) = Ki2{^, t) + Ki^ii, t) , e G r2 , t G (0, T] , (3.44) 

where Ki2{-,t) G I2 H Z and Ki4^{-, t) ^Yt^T] Z are chosen as follows: 

a) The profile Ki2{(,,t) is the unique solution, given by Lemma 13.41 of the equation 

AKi2 + P2{V^^ ■ VF,) + P2(A • VFi) + -dtF^^ = , (3.45) 

ai 

where P2 is the orthogonal projection in Y onto 12- 

b) The profile Ki4{^,t) is the unique solution, given by Lemma [331 of the equation 

AKii + Pi{V^^ ■ VF,) + Pi{Di ■ VFi) + Ci = 0, (3.46) 

where P4 is the orthogonal projection onto Y4, and Ci is as in (13. 7p . 
The main result of this section is: 

Proposition 3.7 Fix ^ < 7 < 1. There exists C > such that, with the above choices of the 
vorticity profiles Fi, Hi, Ki, the residuum of the approximate solution ()3.34p satisfies 

\I^\m<c{^f\-^\^\'l\ (3.47) 
for all^eM."^, allt G (0,r], and alii G {I,-- - ,N]. 

Proof. The proof is a long sequence of rather straightforward verifications, some of which will 
be left to the reader. We first summarize the informations we have on the vorticity profiles Fi, 
Hi, Ki, and on the associated velocity fields u^', v^' . Given any 7 < 1, we claim that there 
exists C > such that 

\Fi{^,t)\ + \VFi{C,t)\ < Ce-^l«l'/^ , \v^^{tt)\ < ^^^^^^,^3^, , (3.48) 

for all ^ G M^, ah t G (0,r], and all i G {1, • • • , A^}. Indeed, we recah that Fi = Fi + F^ , 
where F^ is defined by (j3.39p and Fi by (j3.40p . From (j3.39p and Remark 13.61 we know that 
F^{-,t) G ^2 n Z^. Thus F^ and VF^ satisfy the Gaussian bound in (13:481) . and it follows 
from [18, Proposition B.l] that \v^'{^,t)\ < C(l + |^P)"^/^, see also dSJH). On the other hand, 
using (|3.42p and the explicit expression of the integral kernel of the semigroup S{t) = exp(r£) 
[181 Appendix A], it is straightforward to verify that |Fj(^,t)| < Ce""^'^' Since VS'(t) = 
e'^/^S'(r)V, we also have 

VF,{;t) = -^^*5(log^)(^)"vPoQ.(-,s)ds, 
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from which we deduce that \VFi{^,t)\ < Ce-^l^l'/l Finahy, since Fi{-,t) is a radiahy symmetric 
function with zero average, we have a simple formula for the associated velocity field 



2^ 1^1 J\i='\>m 



l?'l>l€l 

which implies that v^'{£,,t) has also a Gaussian decay as |^| —>■ oo. This proves ()3.48p . 

The corresponding estimates for Hi{^,t) and Ki{^,t) are easier to establish. Since Bi(-,t) £ 
Y^nZ^, it follows from (103]) and Remark ESI that Hi{-,t) G n for all t G (0,r]. Using 
the same arguments as before, we obtain 

\H,iC,t)\ + \VH,iC,t)\ < Ce-^l«l'/4 , < (T^^ • (3-49) 

Similarly, it follows from (lOSl) . ([3M\\ that Ki2{-,t) £ Y2 n and Ki4,{-,t) £ Y^n Z^. We 
conclude that Ki = Ki2 + Ki/^ satisfies 

\m,t)\ + |Vi^,(^,t)| < C7e-^I«IV4 , \v^^{i,t)\ < ^ ^^^^^3^^ . (3.50) 

Next, we make the following observations, which were implicitely used in the definitions of the 
profiles Fi and Ki . Since = + i^^ G + ^2 , it is easy to verify that • VFj G + >2 + . 
Similarly, using the definition (j3.4ip of the vector field Dj, we find that Di ■ VFj G 10 + ^2 + ^4- 
Thus we have the identities 

v^^ ■ VFi = (Po + P2 + Pa){v^' ■ VFi) , A • VFi = {Pq + P2 + Pi){Di ■ V Fi) . (3.51) 

Moreover, it is straightforward to check that 

Po{v^^ ■ VFi) = Pq{v^^ ■ VFn , Po(A • VFi) = Po(A • VFn ■ (3.52) 

Remark that both expressions in (j3.52p appear in the definition (j3.4Up of Fj. 

Now, we replace the definitions (|3.38p ~ (|3.40p and (|3.43p - (j3.46p into the expression (j3.37p of 
the residuum Rf{S^,t). Using in addition (j3.6p . (j3.5ip . (j3.52p . we obtain the simple formula 

ait p(o)/. .X . /i^i 



RKtt) = nritt) + [^) \i^,t) ■VFiiC,t) , (3.53) 
where 

A.«.*)^E?("1^-^)-1^))-f°.(U). (3,M) 

Our goal is to obtain an estimate of the form ()3.47p for Rf{^,t). Since Rl^\^,t) satisfies (13. Sp . 
it is sufficient to bound the second term in the right-hand side of (|3.53p . As in the proof of 
Proposition 13. H we can assume that |^| < d/(2\/i4), because in the converse case the quantity 
|VFi(^,f)| is extremely small due to (j3.48p . Using the decomposition (|3.10p and the definition 
(I3.4ip of the vector field Di , we find 



-.«.')-^i:?(v.(e^).Mf-^))-f°.«-) 
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where 



In view of (|3.1ip . the contributions of V2 are neghgible, and using (|3.48|) we easily obtain 

which is the desired estimate. 

To complete the proof of Proposition 13.71 it remains to verify that 

\Rr\tt)-Ri{tt)\<cQf\~-\^\y^. 

This follows immediately from (j3.36p , (j3.37p if one uses the bounds (j3.48p , (j3.49p , (|3.50p on the 
vorticity profiles Fi, Hi, Ki and the associated velocities. In particular, it is straightforward to 
check that each term in the difference i?^'* — is of the order of [ut/d^)^ for some n > 3/2, 
either due to an explicit prefactor or as a consequence of the polynomial decay of the velocity 
fields v^, v^', or as |^| — > 00. This concludes the proof. □ 

Remark. Instead of (I3.43p . one can define the profile Hi{^,t) as the (unique) solution of the 
elliptic equation 

^^H^-CHi) +AHi + Bi = , 



ai\2 

which is the analog of (I3.39p . In the same spirit, one can replace p.45p by 

-{2K,2 - CKi2) + AKi2 + P2{V^^ ■ VFi) + P2(A • VF^) + -dtF"^ = , 

Oi Oi 

and proceed similarly with (j3.46p . After these modifications, it is easy to verify that the residuum 
satisfies the improved bound 

which is sharper than (j3.47p for small times. This refinement is not needed in the proof of 
Theorem 12. 5|, but it indicates the correct way to proceed if one wants to construct even more 
precise approximations of the A'^-vortex solution. 



4 Control of the remainder 

In the previous section, we constructed an approximate solution t), v'l^^{(,, t) of equation 

(j3.ip . with a very small residuum. We now consider the exact solution u)j(^,t) of (|3.ip given 
by (|2.13p and Lemma [221 and we try control the difference Wi{^,t) — w^'^'^{^,t) for ^ G and 
t € (0, T], in the vanishing viscosity limit. To this end, it is convenient to write 

w.{^,t) = wr''iC,t) + (^) mitt) , v.iC,t) = + (^) Viitt) , (4.1) 

and to study the evolution system satisfied by the remainder Wi{£,,t), Vi{S^,t). 
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Replacing (j4.ip into (j3.ip . and using the definition (|3.4p . we find 

tdtWiiC, t) - {£wi){C, t) + WiiC, t) (4.2) 

+ • V^2;,(C,t) + • V<PP(^,t)) (4.3) 

-E?{"r(«+^-)-"°(^)}-v^-.«.) (-) 

+T.^-^{^+TM^t).vwric,t) (4.5) 

+ E¥^^(^+7=^'*)-^**(^'*)+^*(^'^) = 0' (4-6) 



where Ri{tt) = {vt/(f')-^Rf\i,t). From Proposition [321 we know that 



\Hm<cQ"'\~^\^\''\ (4.7) 

for all i G M^, all t E (0, T], and ah i G {1, . . . , iV}. Also, since !^2Wi{i, t) d^ = 1 by (f2l3]) and 
Lemma [2^21 it is clear that j^2Wi{i, t) d^ = for all t G (0, T]. 

To prove Theorem 12. 5| our strategy is to consider the (unique) solution Wi{£^,t) of sys- 
tem ()4.2p - ()4.6p with zero initial data, and to control it on the time interval (0, T] using an 
energy functional of the form 

N 

where the weight functions Pi{(,,t), i £ {1,...,A^}, will be carefully constructed below. In 
particular, we shall require that Pi{(,,t) > Ce^'^'/^ for some /? > 0. Using (|4.2p - (|4.7p . we shall 
derive a differential inequality for E{t) which will imply that E{t) = 0{vt/<P) as — > 0. This 
will show that 

N 

? X/"^'^V.(e,*)l^de < C(^) , t G (0,r] , (4.9) 



i=l 



for some C > 0, if is sufficiently small. In view of ()4.ip . this estimate is equivalent to (|2.20p . 
which is the desired result. 



4.1 Construction of the weight functions 

Since the construction of the weights Pi{i-,t) is rather delicate, we first explain the main ideas 
in a heuristic way. Ideally, we would like to use for each i G {!,..., A'^} the time-independent 
weight p{^) = e'^l in order to control the remainder Wi{-,t) in the function space Y defined 
in (j3.17p . This is a natural choice for at least two reasons. First, the linear operator £ defined 
in (13. 2p is self-adjoint in Y, and a straightforward calculation (which will be reproduced below) 
shows that 

/ e\^\'/'mi^w.-u>^)d^<-[ e\^\'/'(hvm\'+^S\^^\' + hm\')d(. (4.10) 

Jr2 Jir2 V4 24 2 / 
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Next, using (j3.34p . we observe that the self-interaction terms (j4.3p have the form 



^ t) ■ VwiiC, t) + v,{^, t) ■ V<PP(e, t)) = ^ Kw,{i, t) + regular terms , 

where A is the linear operator defined in (13.190 . Here and below, we call "regular" all terms 
which have a finite limit as — > 0. Since A is skew-symmetric in the space Y by Proposition 13.31 
we see that the singular term {ai/v)K'Wi will not contribute at all to the variation of the energy 
if we use the Gaussian weight e'^' 

Unfortunately, this naive choice is not appropriate to treat the advection terms ()4.4p . which 
describe how the perturbation Wi of the vortex is transported by the velocity field of the 
other vortices. Indeed, using again ()3.34p . we can write (14. 4p as 

y ^LG/^^£i0\ _^G/£^\1 + regular terms . (4.11) 

Since is given by (|1.7p . and since > d > when i / j, it is easy to verify (as in 

the proof of Proposition 13. ip that the first term in (j4.1ip has a finite limit as ^ 0, provided 
|.^| ^ dj^fvi. Although this is not obvious a priori, we shall see below that the same term is 
also harmless if \(\ ^ D/\/ljt, where 

D = max max IzJt) — Zi(t)\ < oo . (4.12) 
te[o,T] i^j ' ^ ' J"- 

In the intermediate region, however, the first term in (j4.1ip can be of size Odal/i/), and there 
is no hope to obtain a better bound. But we should keep in mind that the whole term (|4.4p 
describes the advection of the perturbation Wi by a divergence-free velocity field, and therefore 
does not contribute to the variation of the energy in the regions where the weight function 
is constant. The idea is thus to modify the Gaussian weight to obtain a large plateau in the 
intermediate region where the advection term (|4.4p is singular. 

Our second try is therefore a time-dependent weight of the form 



el«l/4 if \^\<p{t), 
if < < 

,\e/i^K^) if 1^1 > Kp{t) , 



ePWV4 if p(^t) < 1^1 < Kp{t) , (4.13) 



where 



d , 4L> 

P^*) = TTT^ ' and K = — 



By construction, the function ^ 1-^ p{^, t) coincides with e'^' in a large disk near the origin, is 
identically constant in the intermediate region where the advection terms (|4.4p are dangerous, 
and becomes Gaussian again when |^| is very large. Moreover, this weight is continuous, radially 
symmetric, and satisfies el€l'/(^^') < p{^,t) < el«l'/^ for all ^ G and all t G (0,r]. Of 
course, the operator C is no longer self-adjoint in the function space defined by the modified 
weight p{(,,t), but an estimate of the form (j4.1Up can nevertheless be established by a direct 
calculation. Similarly, the operator A is no longer skew-symmetric, but we shall see below that 
it remains approximately skew-symmetric with the modified weight, and this will be sufficient 
to treat the self-interaction terms (14.31). 
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We now consider the contributions of the advection terms (j4.4p to the variation of the energy 
(I4.8p . when pi{^, t) = p{^, t) for alH G {1, . . . , N}. Integrating by parts, we obtain 



If pit) < \(,\ < Kp{t), then Vp(^, t) = by construction, hence it remains to consider the regions 
where |^| < p{t) or |^| > Kp{t). Unfortunately, although the integrand in (|4.14|) is regular when 
ICI ^ p{t)i the contributions from that region are still difficult to control if t is large. To see 
this, we first replace v^"^ by v'^ as in (|4.11|) . because the difference v^"^ — v'^ is negligible at 
this level of analysis. Assuming that |^| < we find as in (j3.54p 

where is defined in ^Al\i . As Vpi(,,t) = for |^| < p{t), we conclude that the 

main contribution to (j4.14p has the form 



4tt (P 



I fetetji^^a -'i^^^WjK^^.w (4.15) 



'l«l<P(t) 

and is therefore bounded by 

/ |^|^el«lV4|^^(e,t)Pde. (4.16) 

If t is small with respect to the turnover time Tq defined in (j2.10p . then \ai\t <C and the 
quantity (j4.16p can be controlled by the negative terms originating from the diffusion operator 
see ()4.10p . In that case, one can show that the expression (|4.14p is harmless also in the outer 
region where |^| > Kp{t), and it is not difficult to verify that the linear terms (14. 5[) and the 
nonlinear terms (|4.6p can be controlled in a similar way. Thus, if T ^ Tq, it is possible to carry 
out the whole proof of Theorem 12.51 using the energy functional (j4.8p with Pi{^,t) = p{^,t) for 
alii G {1,...,A^}. 

The main difficulty is of course to get rid of the condition T <^ Tq, which is obviously too 
restrictive. We follow here the same strategy as in the construction of the approximate solution 
w^^^{S,,t) in Sections 13.21 and 13.31 So far, we have used radially symmetric weights to eliminate, 
or at least to minimize, the influence of the singular self-interaction terms (14. 3p . The idea is now 
to add small, nonsymmetric corrections of size 0{iyt/d?) which, by interacting with the singular 
expression (14. 3|) . will produce counter-terms of size 0(1) that will exactly compensate for ()4.15p . 
Unfortunately, the construction of these corrections is quite technical, and requires a non-trivial 
modification of the underlying radially symmetric weight. 

We now give the precise definition of the weights Pi{^,t) that will be used in the definition 
(j4.8p of the energy. We start from the radially symmetric weight 



el5lV4(^l_^(|^|2_a(i)2)j if < < a(t)2+l 
e'^W'/^ if a(t)2+l < |^|2 < 6(t)2 , (4.17) 

e/5|«l/4 if |^P>6(t)2 
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where 



a(t) = ao(-) , K*) = ^o(-) 



and ^ = -0 < 1 
bo 



(4.18) 



Here oq <C 1 and 6o ^ 1 are positive constants which will be chosen later. In ()4.17p . we use a 
cut-off function : (— cxd, 1] E which satisfies V'(y) = for y < —1, ^(—1) = 1) = 0; 
< ■(/''(y) < 3(1 - ^{y)) for \y\ < 1, = 1 - e~^/'*, and ^'(1) = 36"^/"^. For definiteness, we 
can take 

V'(y) = C2(y + i)'-C3(2/ + i)' for |y|<i, 

where C2 = | ~ 'E^~^^^ ~ 0.068 and Cs = 3 ~ ~ 0.0066. It follows from these definitions 

that the function ^ 1-^ Po{S,yt) is piecewise smooth and nondecreasing along rays. Moreover, we 
have 



' |eieiV4(^l_^(|^|2_«(i)2)^ if < < a(t)2+l , 
ff a(t)2+l < < 

if |eP>6(t)2, 



Vpo(^,t) 

where V'(y) = V'(y) + The graph of the function po{C,t) is depicted in Fig. 1. 



(4.19) 



Po(^,i) 



=/3|«l/4 



a(t) 



Kt) 



Fig. 1: The radially symmetric, time-dependent weight po{^,t) is represented as a function of |^|. The gradient 
of po(C)*) has a jump discontinuity at |5| = b{t), but is Lipschitz continuous near |^| = a{t). 



Next, for each i G {1, . . . , N}, we define 
where the correction qi{^,t) vanishes identically when > a(t)^ + 1 and satisfies 



(4.20) 



1 v-^ a,- d'^ 



^^(0-vg.(e,t) = - E- 



vr ^ a, |2:ij(t)|' 



(4.21) 



for < a(t)^ + 1. Proceeding as in Section [3l it is easy to obtain the explicit expression 



1 e-Vpo(e,i) v-oj- d 



2 i_e-ICP/4 Z^a. 1^ (t)|4 



(^z,,(^))^-(^z,,(^)^)^ , (4.22) 
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which is in fact vahd for |^| < b{t), because 

Vpo(^,i) vanishes when a{tf + 1 < < b{tf. It 

fohows from KWl that 

{l + \e)m,t)\ + m'^Q^im + \t^tq^{tt)\ < c\c\\i + m'poic,t) , (4.23) 

whenever < a{tf + 1. Using K20\\ and the definition (^TJEt} of a{t), we deduce from K23\\ 
that \pi{(,,t) — Po{S,jt)\ < CaQPo(C)^) for some C > 0. Thus, if we choose the constant ao > 
sufficiently small, we see that the weight Pi{^,t) is a very small perturbation of po(^)0 foi" ^-H 
^ G R2 g^ij f g jj^ particular, for i = 0, . . . , A^, we have the uniform bounds 

i e^l«l/^ < p^it t) < 2 el«l'/^ , C&R\ i G (0, T] . (4.24) 
4.2 Energy estimates 

Now that we have defined appropriate weights ,PN{^,t), it remains to control the 

evolution of the energy functional E{t) introduced in (j4.8p . Here and in what follows, we always 
assume that the parameters in ()4.18p satisfy oq <C 1, 6o ^ 1, and that the viscosity > is 
small enough so that uT <^ d?. 

Proposition 4.1 There exist positive constants eq, ei and n\, K2, K3, depending only on the ra- 
tio T/Tq, such that, if = 6q"^ = eo and vT/d'^ < ei, and if w = {wi, . . . ,10^) G C^{{0,T],Y^) 
is any solution of system (|4.2p - ()4.6p . then the energy E{t) defined in (14. 8p satisfies 

tE'{t) < -KiE{t) + (^(t) + Eitf'^ + (^^) , < t < r . (4.25) 
Proof. Differentiating (j4.8p with respect to time and using (|4.2p - (|4.6p . we find 

N 6 

tE'{t) = / ,{-t5tPi(e,t)|iI;i(e,t)P+Pi(e,t)u'i(^,t)t5tu;i(e,t)}de = J]^fc(t) , (4.26) 
i=i -^^^ k=i 

where 

Slit) = Y] / |7T*5tPi(e,t)i^i^i(e,i)P+Pi(e,t)tSi(c,i)(>c^s^(e,t)-tii(e,i))|de , 

iM2 L 2 J 

N 

-^2(0 = -E / p.(e,i)*«(^,o-(<""(e,t)-VTi),(c,t) + u,(e,t)-v<pp(e,t))dc, 

^3(0 = - E E ? {-r(f + ^ • «) - } • " • 

^4(i) = - E / t)M^, t) E - (e + ^ , t) • vu;r(^, o d^ , 

AT » ^ / (f) 

= - E / Pi(e,t)«)i(c,oE^^^(^+%^'0-^*^(^'*)<i^' 

N 

£6{t) = "E / PiiC,t)Wi{^,t)R,{^,t)dC . 



28 



The general strategy is to control the contributions of £2(^)1 ■■ ■ j^eit) using the negative terms 
contained in £i{t). We shall proceed in six steps, each one being devoted to the detailed analysis 
of a specific term. Very often, we will have to consider separately the three cases < a(t)^ + 1, 
a(i)^ + 1 < 1^1^ < b{t)^, and |^| > b{t), see (I4.17p . For simplicity, these three domains in will 
be denoted by "region I", "region II", and "region III", respectively. 

Step 1: Diffusive terms 

Our goal is to show that there exists a constant k > such that 



N 



(4.27) 



i=l 



where 



PoiC,t)\Vwiic,t)\^dC , Ji{t)= / Mi,t)x{i,t)\wi{i,t)\''di 



it) = [ Po{^,t)\m{tt)\' 



(4.28) 



Here t) is the continuous, radially symmetric function defined by 

r \e if < lei < ait) , 

Xi^,t) = \ a{tf if a(t) < lei < 6(t) , 
I m\ if lel > hit) . 

To prove (j4.27p . we start from the identity 

/ piWi{Cwi -Wi)d^ = - I (pi\Vwif + Wi{Vpi ■ Vwi) + ^(^ • Vpi)\wi\'^ + l:Pi\wi\'^] d^ 
Jr2 Jr2\ 4 2 / 

which is easily obtained using (13.21) and integrating by parts. Since 



Wi{Vpi ■ Vwi) d^ 
we see that £i{t) < £i{t), where 

N 



< 7 / Pi|V?ii|2de + 



1 f \Vp. 



*l l„r. |2 



R2 Pi 



Wi\ d^ , 



= - E / (7K(e,t)|vuii(e,t)p+K(e,t)x(e,t)i«^i(e,t)i' + jK(e,i)i^^i(e,t)i')de 

^7r2V4 2 / 



and 



m,t) 



^■vpiitt) ivp,(e,t)|2 tdtPi{^,t) 



4p*(e,t) 3pi(e,t)2 2p^{tt) 

We already observed that Pi{^,t) > ^Po{£,,t) if uq is sufficiently small. So, to prove (j4.27p . it 
remains to show that x(C) t) > Cxi^,, t) for some C > 0. This is easily verified in regions II and 
III, because Pi{^,t) = po{C,t) in these regions. From (j4.17p . we find by a direct calculation 



xiC,t) 



16 

16' ' 48 



if a(t)2 + 1 < < 6(t)2 , 
if \^\>b{t). 
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In region I, we first compute the contributions of the radiahy symmetric weight Po{^,t) to the 
function xiCii)- Using (I4.19p . we obtain 



8 l-V(|CP-a2) 12 Vl-V(|eP-a2)/ 4 1 - V^d^P-a^) ' 

In particular, we see that xo(C,0 = |C|V(24) if |CP < a{tf - 1. If y = |Cp - a{tf G [-1,1], we 
use the elementary bounds 

0<i4M<:, and i(i^iM)±MM > ^ > . 

l-V(y) l-V'(y) 

which follow from the definition of ip, and we deduce that Xo(C,i) ^ 5min{|^p, a(t)2} for some 
(5 > 0. Summarizing, we have shown that there exists Co > such that xo(C, ^) ^ C'ol^P in region 
I. On the other hand, using ()4.23p . it is straightforward to verify that ^ Xo(C, "C'iOqI^P 
when |^|2 < a(t)2 + 1. Thus, if the constant oq in ()4.18p is sufficiently small, there exists C2 > 
such that x(C, t) > C2X(?, for ah ^ e and all t E (0, T]. This concludes the proof of (lOTl) . 

Step 2: Self-interaction terms 

We next consider the term <f2(i) in ()4.26p . To simplify the notations, we rewrite (j3.34p in the 
form 



wnc,t) = + [^)m,t) , = t^^io + [^y-i^^t) > (4-29) 

where = Fi{^,t) + (i^t/ff)^/^ Hi{^,t) + {i^t/(f)Ki{C,t), and v-^' is the velocity field ob- 

tained from J^i via the Biot-Savart law ()1.3p . We thus have <?2(i) = f^i(t) + • • • + f^4(i), where 



N 
N 



i=l 
N 



To prove that has a finite limit as — > 0, we integrate by parts and use the fact that 

v'^ • Vpo = because the weight po is radially symmetric. In view of ()4.20p . we find 



i=l -"^ i=l 



This term cannot be controlled by (j4.27p . unless t <^ Tq, but it will exactly compensate for 
another term coming from £3{t). As was already explained, this is precisely the reason why 
the correction {h't/(P)qi{(,,t) was added to the weight Po{(,,t). To treat ^2, we use the fact 
that the linear operator Wi Vi ■ VG (where Vi is obtained from Wi via the Biot-Savart law) is 
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skew-symmetric in the function space Y defined by (j3.17p . see [19]. We thus have 



\^2{t)\ 



N 

~1 Jr^^ ' 



i=l 
N 



<Y, — I e\^\''^\w,{i,t)\\vi{i,t)\\VG{i)\di 

< C^:^ / e/^l«l/V.(e,t)||^.(C,i)||^|e-^l«l/«de . 

In the last hne, we have used the definition p.7p of G. To estimate the last integral, we apply 
the trilinear Holder inequality with exponents 2, 4, 4. We have ||e^l^l/^tDj||/,2 < CK^J"^ by (j4.24p . 
()4.28p . and ||'Uj||j;,4 < C||7l;i||j;^4/3 < G/C/ (see e.g. [S]). Moreover, by direct calculation, 

N 1/4 p-/3a{t)/8 
l?P>a(t)2-l ^ (/3a(t))V4 

provided that /3a(t) > 1. Thus, if the constant ei in Proposition 14.11 is sufficiently small (de- 
pending on eo and T/Tq), we find 

AT I I N 

mt)\ < cE^K:,(t)a(t)3/2e-^"W/« < eY,m) ■ 
1=1 1=1 

Here and in what follows, e denotes a positive constant which can be made arbitrarily small by 
an appropriate choice of eo and ei. Using similar estimates, it is also easy to bound the regular 
terms J^s and 1^4. We find 

mt)\ < y^/o(c,t)i«)^(e,i)i(k^"(e,t)iivu;i(e,t)i + i^.(e,t)iiv.F,(e,t)i)d^ 

TV I I N N 

i=l i=l " i=l 

Step 3: Advection terms 

Using (j4.29p and integrating by parts, we write <?3(t) = '^lif) + ^2{t) + ^3(i)) where 



- E/j-«.*)Pv.„,u).i;g{.<^(«.af ) - .°(2f )} . 

7^1 ^2d2 V ^ y 
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By construction, all three integrands vanish for a(t)^ + 1 < |Cp < ^(0^) so we need only consider 
regions I and III. The most important term is the contribution of region I to ^i{t), which we 
denote by "^{{t). Using (j3.54p together with the identity 

v'^io ■ vgi(e, t) = -vpoic, t) ■ A(e, t) , lep < a(t)2 + 1 , 

which follows immediately from the definitions (j3.4ip . (j4.2ip . we find 

N 

^f(t) = + / |u;,(e,t)|2vpo(e,t)-A,(e,t)de. 

^ ^^|CP<a(t)2+l 

Thus, using the estimates on Aj(^,t) which follow from (j3.55p . (j3.56p . we obtain 

i=l 4=1 

On the other hand, the contribution of region III to ^'i(t) can be estimated as follows 

N 



bo f-^J\£\>b(t) 



=1 J\<i\>h(t) 

if the parameter 6o > is chosen sufficiently large. Similarly, using ()4.23p . we can bound 
in the following way 

AT 



if ao > is sufficiently small. Finally, to bound ^'3(t), we recall that \v^'{tt)\ < C(l + |Cp)"^/^ 
see dSaSD, dSaiD, (|330l) . Since |^ + ^| > ^ in regions I and III when i / j, we find 

Step 4: Cross-interaction terms 

To bound i?4(t), we need a good estimate on the velocity field Vi{S,, t). We claim that there exists 
a constant C > such that 



ii(i + iei')i)i(e,t)iiioo < c(j,(t) + /c,(t)j , ie{i,...,iv} . (4.30) 

Indeed, since J]g2?^j(^, i) d^ = for all t G (0,T], it follows from [18, Proposition B.l] that 
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As H^M.^) ^ L^{m^) and (1 + < Cpo{^,t), we have ||(1 + \^\^)wi\\L4 < C{ll^^ + IC]^^). In 
addition, using Holder's inequality, we find ||(1 + |'^P)ii'i||i4/3 < C/C|^^. This proves (|4.30p . 

The main contribution to £4{t) comes from region I. Since \^ + > 'j— when i j and 
I^P < a{tf + 1, we obtain, using (|T30]) . 

\£l{t)\ < C5]5]^fc(t) + /C,(f))'^' / Poimmi^M^w^^'imdC 

I I TV 1/2 ^ ^ + ^ 

j=l i=\ i=\ ^ i=\ 

In regions II and HI, the quantity | Vw^^''(^, t)| is bounded by Ce~^^^^^^ for any 7 < 1. Choosing 
7 > 1/2 and proceeding as in the second step, we easily find 

AT I I N 

\el\t)\^\£i'\t)\ < c7^^^x;.(t)V%.(t)i/2e-(7-|wt)V4 < e^icm ■ 

Step 5: Nonlinear terms 

Instead of (j4.30p . we use here the simpler inequality 

2 ^ 1/2 

WViWl^ < C\\Wi\\L4\\Wi\\^4/3 < C[Ii+ICij /c/ 



which follows from |18l Lemma 2.1]. Applying Holder's inequality, we thus obtain 

N N 

i=i ° i=i 

Step 6: Remainder terms 

Finally, using Holder's inequality and estimate (j4.7p . we find 

4=1 i=l 

Collecting the estimates established in Steps 1-6 and using the fact that E{t) « i J2iLi 
we easily obtain (j4.25p provided that e > is sufficiently small. As was already explained, this 
last condition is easy to fulfill if we choose the constants eo, ei appropriately. This concludes 
the proof of Proposition 14.11 □ 

4.3 End of the proof of Theorem 12.51 

It is now quite easy to conclude the proof of Theorem 12.51 If the solution uj'^{x,t) of ()1.2p with 
initial data (j2.ip is decomposed as in (j2.6p . we know from Lemma [2.21 that the rescaled vorticity 
profiles Wi{^,t) = defined by (f233|) satisfy Wi{-,t) G Y for all t G (0,r], see (USD- 

Standard parabolic estimates then imply that w = {wi, . . . , W]\f) E C^{{0, T], 1"^) is a solution 
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of system (j3.ip . where the vortex positions Zi{t) = z^{t) are given by (j3.4p . Moreover, we know 
from [151 Proposition 4.5] that Wi{-, i) — > G in y as t — > 0, for ah i G {1, . . . , N}. In fact, in 
Section 6.3 of [15], this convergence is estabhshed in a polynomiahy weighted space only, but 
the proof also works (and is in fact simpler) in the Gaussian space Y. Using the approximate 
solution w^^'^{(,,t) of (j3.ip constructed in Section [3l we can even obtain the following improved 
estimate for short times: 

Lemma 4.2 For any fixed u > 0, one has \\wi{-,t) - Wi^^{-,t)\\Y = 0{t^/^) ast ^ 0, for all 
ie{l,...,N}. 

Proof. Let Wi{^,t) = Wi{S,,t) — w^'^'^{£^,t). Then Wi satisfies a system which is similar to (I4.2p - 
(|4.6p . except that the first line (j4.2p is replaced by tdtWi{£,,t) — {C'Wi){S,,t), and in the last line 
(14.61) the left-hand side becomes 

N .X 

^^U^+'-im,t). VMC, t) + Rf (^, t) , (4.31) 

/ON 

where Ri {£,,t) satisfies p.47p . Note that the nonlinear terms in (j4.3ip are singular as ^ 0, 
but this is not a problem here because > is fixed. To estimate Wi{-,t) in the space Y, we 
use the energy functional 

Repeating the proof of Proposition 14. 1^ with substantial simplifications, we obtain for small 
times a differential inequality of the form 

tE'{t) < -viEit) + V2E{tf + m (^) ' , (4.32) 

where the positive constants 77i,??2,f/3 may depend on v. In the derivation of (j4.32p . the only 
new ingredient is the estimate 

/ e\^\'/^mCmd^ < -k[ e\i\'/^(\Vm\^ + \C\^\m\^ + \m\AdC, 

which holds for some k > because the self-adjoint operator £ is strictly negative in the subspace 
of functions with zero mean, see [IB]. Since we already know that £'(t) — > as t — > 0, inequality 
([02|) imphes that E'{t) < Ct^ for t > sufficiently small, hence E{t) = 0{t^) as t ^ 0. This 
is the desired result. □ 

We now consider the energy functional E{t) defined by (|4.8p . Since Wi{S,j t) = {ut/ d?')wi{S,, t), 
and since the weights Pi{i,t) satisfy ()4.24p . it follows from Lemma 14.21 that £'(t) — > as t — > 0, 
for any fixed u > {). As long as E[t) < 1, we have by Proposition 14.11 

tE'it) < 2k2 ^^E{t) + 

hence 

Eit) < e^-^^/^o (^) ds < e^-VTo . (4.33) 

If we choose v > Q sufficiently small so that 

we see that ([^1) holds for ah t € (0, T]. Since pi(C, t) > Ce^l^l/^ for all ^ G ^^^^ ^ ^ (g, T], 
we obtain ()4.9p . and (j2.20p follows. The proof of Theorem 12.51 is thus complete. □ 
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5 Appendix 

Proof of Lemma 12.21 Let U : M? x (0, oo) ^ be a smooth, divergence-free vector field, 
and fix > 0. We assume that 

-SUp( {ut)^/^U{;t)\\L^ + m;t)\\LA = Cq < OO , (5.1) 

where = d\\J2 — d^^x- Then any solution of the linear equation dt^o + {U ■ V)ti; = z^Ao; can be 
represented as 

w(x,t) = / rf}(x,t;y,s)u;(y,s)dy , x G , t>s>0, (5.2) 

where the fundamental solution r|^(x, t; y, s) has the following properties: 

1. For any P £ (0, 1), there exists Ci = Ci(/?, Co) > such that 

< rU^. t: y, s) < exp(-3 , (5^3) 

for all X, y G and all t > s > 0. This very precise upper bound is due to E. Carlen and 
M. Loss [1]. 

2. There exists 7 = 7(Co) > and, for any 5 > 0, there exists C2 = 02(6, Cq) > such that 

\T'{j{x,t;y,s)-T'{j{x',t';y',s')\ < C2[\x - x'\^ + \t - t'\^/^ + \y - y'\^ + \s - s'l^y^^) , (5.4) 

whenever t — s > 6 and t' — s' > 6. This Holder continuity property, which is due to H. Osada 
[15], implies in particular that T'^{x, t; y, s) can be continuously extended up to s = 0, and that 
this extension (which is still denoted by T^) satisfies ()5.3p and (j5.4p with s = 0. 

3. For all x,y gE? and alH > s > 0, we have 

/ T'{j{x,t]y,s)dx = 1 , and / T'{j{x,t;y, s) dy = 1 . (5.5) 

Note that the first equality uses the fact that U is divergence-free. 

We now consider the particular case where io{x,t) = uj'^{x,t) and U{x,t) = u'^{x,t). Then 
dtLV + {U ■ V)w = I'Auj by construction, and the results established in [21] show that assumption 
(15. ip is satisfied with Co = C\a\/u, where C > is a universal constant and \a\ = |ai| + • • • + 
|QAr|. Using the Holder continuity (j5.4p and the fact that uj'^{-,t) ^ /i as t — > 0, we can take the 
limit s ^ in the representation (j5.2p and obtain, for any > 0, the following expression 

N 

uj''{x,t) = / r^.(x,t;y,0)d/Uj; = air^,.(x, t; Xj, 0) . (5.6) 

Setting wj" (x, t) t;xj,0), we obtain the desired decomposition ()2.6p . and the various 

properties of u;f(x,i) follow directly from (15. 3p and (15. 5p . □ 

Proof of Lemma 12.31 For simplicity, we set z'^{t) = {z'({t), . . . , z'^{t)) G (M^)^ and we rewrite 
system (i2TTD as z^it) = F {z" {t) , ut) , where F : (R^)^ x (0,oo) (M^)^ is defined by 

^^(-'^) = E^-^(^)' ^^{l,■■■,N}. (5.7) 
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For any 5 > 0, we denote = {z G (M^)^ ; \zi — Zj\ > S for all i 7^ j}- Then F extends to a 
smooth map from no X [0,00) to (R^)^ , with 

F'^i^^^)=T.?¥^^^ zG{l,...,iV}, (5.8) 



'J I 



This remark already implies that system (j2.1ip is locally well-posed for all initial data in f^o- 
Global well-posedness easily follows, because as soon as r/ is bounded away from zero, the vector 
field z I— > F{z, rf) is smooth and uniformly bounded. 

We now compare the solutions of (j2.1ip in Qq with those of system (j2.3p . which can be 
written as z{t) = F(z{t),0). z £ for some 6 > 0, then using (j5.7p . ()5.8p and the definition 
(|1.7p of v'-^, we easily find 



\F,iz,rj)-F,iz,0)\ < J^M^^e-l^-^^l'/(4^) < 



27r Izj — Zi\ 27:6 



|a| _a2 
■ e 



52 /(4»;) 



for any 77 > 0. Similarly, if z, z £ Qg, then 

\Fi{z,0) - Fi{z,0)\ < ^ max{|zi - zi|, . . . , |zAr - zn\} < - z\\ , 

where \\z — z\\ = max{|zi — zi\, . . . , \z]y — zj\f\}. 

Assume now that z G C0([0, T], (M^)^) is a solution of ([23]) satisfying i^M) for some d > 0, 
and take 5 € (0, d). For any > 0, let z'^{t) denote the unique solution of (j2.1ip with initial 
data z'^{0) = z{0). As long as z'^{t) stays in 0^, we have 



P"(*) - MM < \lF(z-'{t),i,t) - F(z''(*),0)|| + ||F(=-((),0) - F(=((),0)|| 

hence 



2vr5 Jq 2 

If > is sufficiently small, this implies that z'^{t) G for all t G [0, T], hence (j5.9p holds 
for t G [0,r]. Choosing for instance 6 = d\J 4/5, we obtain (I2.12p with K\ = exp(CT/To). For 
larger values of the solution z^(t) may leave fi^, but in that case the bound (j2.12p still holds 
if we take the constant K\ large enough. □ 

Proof of Lemma 13. 2L Let r = |^|/|7?| < 1, and ip = 9 — (p- We have 

|^ + ry|2 = |77|2(l +2rcos(V') +r2) = |r?|2|l + z|2 ^ 

where z = re*^ G C. Now 

^ ''1 - z + z'^ - z^ + . . .^(l - z + z'^ - z^ 



\l + z\^ 

= 1- {z + z) + {z"^ + zz + z^) - {z^ + z^z + zz^ + f") + 



But, for each n G N, 



z" + z^'^z H h zz""-^ + = = r 



z — z sin(V') 
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if sin('0) 7^ 0. Thus 



|ryp |?7p(|l + zp ) I^P 1^1" sm{ip) ^ 



\^ + v\ 



Multiplying the first and the last member of (15.101) by ^ • r/"*" = \^\\r]\ sin(^), we obtain ()3.12p . 
This concludes the proof. □ 
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